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We present a new QCD description of the ATLAS jet veto measurement, using the Banfi- 
Marchesini-Smye equation to constrain the inter-jet QCD radiation. This equation resums emissions 
of soft gluons at large angles, at leading-logarithmic accuracy, and accounts for both the so-called 
Sudakov and non-global logarithms. We show that this approach is able to reproduce, with no fitting 
parameters, the fraction of high-pr forward/backward di-jet events which do not contain additional 
hard emissions in the inter-jet rapidity range. We also compute the gap fraction in fixed-order 
perturbation theory to 0(a 2 s ) and show that the perturbative series is unstable at large rapidity 
intervals. 



I. INTRODUCTION 



Recently the ATLAS collaboration measured, in proton-proton collisions at the LHC, the fraction of di-jet events 
that do not contain additional hard radiation in the inter-jet rapidity range pQ. The original goal of this measurement 
was to look for BFKL-type effects, as was previously done at the Tevatron with the so-called 'jet-gap-jet' observable 
PHS]- However the use of a veto scale E out ^> Aq CD by ATLAS, instead of a true rapidity gap void of any hadronic 
activity, drastically reduces the sensitivity to BFKL physics. Rather, it turns out that this 'jet-veto' measurement is 
sensitive to the physics of inter-jet energy flow, and in the limit pr 3> -Eout, to the resummation of soft large-angle 
gluon emissions in perturbative QCD. 

In that regime, the picture that emerges from the comparison [T] between the ATLAS data and Monte-Carlo 
predictions is far from clear: while HERWIG [5] and PYTHIA [7] tend to be in reasonable agreement with the 
measurement, fixed-order [HI E] calculations matched with parton shower are quite below the data. Furthermore, 
adding BFKL logarithms, as done e.g. in HEJ [ID] , to take into account the fact that the rapidity difference between 
the two primary jets, Ay, can be large does not help. A possible origin of these shortcomings may be the insufficiency 
of the angular-ordered parton shower to faithfully capture the physics of energy flow. Indeed, the soft gluons which 
constitute inter-jet radiation are not ordered in angle but rather in px, which makes the relevant resummation single- 
logarithmic (a s \np T /E out ) n . 

The main goal of this paper is to provide a QCD-based computation of the jet veto cross-section obtained by 
resumming these logarithmically-enhanced (a s Inpx/ E out ) n termfj^] There are actually two types of such logarithms — 
the Sudakov logs and the non-global logs. The former is more well-known and can be resummed by using the soft 
anomalous dimension technique [111 112] . On the other hand, the non-global logs arise from soft emissions from the 
secondary gluons (not from the primary hard partons), and can only be resummed in the large- N c limit [13] where 
the successive emission of soft gluons may be viewed as the splitting of color dipolesQ While HERWIG and PYTHIA 
can partly account for the effect of these logarithms due to some overlap in the phase space [15] . in principle they are 
not optimized for observables such as the gap fraction. 

In this paper we investigate whether the ATLAS data can be described by a perturbative framework which incorpo- 
rates the relevant single-logarithms to all orders. For the resummation of the non-global logs, we follow the approach 
of Banfi, Marchesini and Smye (BMS) [16] who reduced the problem to solving a nonlinear integro-differential equa- 
tion. In a previous paper [T7], two of us numerically solved this equation and estimated the survival probability 
of the BFKL-induced rapidity gap. Here we consider the one-gluon exchange (octet) channel and compute the gap 
fraction in proton-proton collisions to be directly compared with the ATLAS data. An approach similar in spirit was 
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The tact that ATLAS quotes very small non-perturbative effects (smaller than 2%) on their measurement of the gap fraction means 
that we may hope to achieve a good description of the jet veto cross-sections based purely on perturbative QCD. 
2 This process is mathematically identical to the gluon splitting in the BFKL evolution [14] . 
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taken in [TH] where the authors did not use the large- N c approximation and fully included the Sudakov logs. The 
non-global logs, however, were included only by introducing a if-factor. For other recent discussions of energy flow 
and the ATLAS data, see Q35H22]. 

The plan of the paper is as follows. In Section II, we first compute the gap fraction in fixed-order perturbation 
theory to 0{a 2 s ) and point out a problem which arises when the rapidity gap Ay becomes large. In Section III, we 
describe the BMS approach which incorporates resummation and compute the gap fraction in this framework. The 
results are then compared with the ATLAS data. Finally Section IV is devoted to conclusion. 



II. FIXED-ORDER COMPUTATION 



The observable of interest is the gapped fraction in proton-proton collisions measured by the ATLAS collaboration 
at y/s = 7 TeV. It is defined by 

K(Ay p ) - I (1) 

' dAy cPpt I dA yd 2 px ' 

where er mcl is the inclusive cross section of di-jet events and Ay is the rapidity difference of the two jets which have 
mean transverse momentum px — (pt,i +Pt,2)/2 > 50 GeV and rapidity \y^\ < 4.4. Jets are reconstructed using the 
anti-fcj algorithm |23j with the radius parameter R — 0.6. In defining the di-jet system in each event, ATLAS used 
two different selection criteria: The highest-pr jet pair and the most forward/ backward jet pair. er vcto is the gapped 
cross section in which a veto is applied to the di-jet cross section requiring that no jet with px above E out — 20 GeV 
is observed in the rapidity interval between the two jets. 

As explained in the introduction, one expects in the perturbative computation of the ratio ([!]) large logarithms 
of the form (ct s lnpx/E out ) n to become important. In particular, this includes the Sudakov-type logarithms 
(a s Ay hipx/E ou t) n which are additionally enhanced by factors of Ay when Ay is large. A fixed-order computa- 
tion should then break down and the large logarithms should be resummed in order to obtain a correct description of 
the vetoed di-jet cross section. However, before coming to that resummation, it remains interesting to check around 
which value of Ay the perturbative scries becomes unstable. Therefore, in this section we first perform a fixed-order 
study of the veto fraction. 

Since the leading-order (0(a 2 )) vetoed and inclusive cross-sections are equal, the first non-trivial order for TZ is 
when the cross-sections are computed at NLO. At first sight, the fact that the inclusive jet cross-section is not known 
at NNLO, because of its 2-loop pure- virtual contribution, seems to indicate that the story ends here and TZ can only 
be computed at 0{a s ). But, since this 2-loop pure- virtual contribution to the jet cross-section is actually the same in 
the vetoed and inclusive cases, one can show that, following the same arguments as in [2"3] , this unknown piece does 
not contribute to the 0{a 2 s ) term in the series expansion of TZ. An explicit series expansion shows that one can write 
TZ a£] 



TZ(Ay,p T ) = 1 - V:: dA !( LO a s 



\^ dp T dAy J 



LO 



\dp T dAy J LO \dp T dAy) NLQ \ 



d PT dAy J NLQ 

Jp T dAyJ LO \dp T dAy) LO 



< (2) 



where the a s dependence has been explicitly factored out and er novcto = cr mcl — er voto . 

All the differential cross-sections in (2) can be computed explicitly using NLOJct++ [25 . Following the ATLAS 
analysis [T] as described above and focusing, for definiteness, on the case where the di-jet system is made of the 
most forward and most backward jets, we simulated events with NLOJct++ v4.1.2 using the anti-/ct algorithm as 
implemented in Fast Jet [26] . In the perturbative computation we have set both the renormalization and factorization 
scalesj^to p T and we have used the MRST2002 PDF set at NLO [27] , 

The result of this analysis is shown on Fig. [TJwhere TZ is plotted as a function of the rapidity interval Ay for different 
bins in px- The dashed (red) curves are obtained by truncating ^ at order a s while the solid (blue) curves also 
keep the 0(a 2 ) terms. We immediately notice that the LO result becomes negative when Ay becomes large while 
the NLO result becomes larger than 1. These unphysical results are actually not surprising and stem from the fact 
that at large rapidity interval Ay one expects a Sudakov-type behavior of the form TZ cx exp(— Ca s Ay \npx) which, 



3 An alternative and maybe more direct approach is to write 1Z = 1 — (j Iloveto /<j mcl and to realize that taking both cross-sections at NLO, 
i.e. (j Ilovcto a t 0(«f ) and cr mcl at C(«f ), would formally give a description of TZ at 0(a 2 3 ), while ct^J^q" does not depend on the 2-loop 
cross-section with only 2 partons in the final-state. The expression in |2| is no more than an explicit series expansion of that statement. 

4 Although this would naturally apply for a di-jet system based on the two hardest jets in the event, px is the natural scale choice also for 
this forward-backward configuration which is expected to be dominated by a t-channel exchange. Note also that we could have studied 
the scale uncertainty and compared it to the differences between the LO and NLO curves in Fig. [Tl 
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Figure 1. LO (0(a s )) and NLO O(a^) predictions for the fraction of di-jet events with a jet veto. 



truncated at fixed order, would give 1Z oc 1 — Ca s Aylnpx + (Ca s AylnpT) 2 /2 + • • • . More interestingly, one sees 
explicitly from Fig. [I] that the perturbative series becomes completely unstable for Ay > 2, with a weak dependence 
on px- The main message of this study is that one should thus expect some form of resummation to be necessary for 
Ay > 2 which is practically the whole interesting region of the physics of rapidity gap. 



III. RESUMMED COMPUTATION 



In this section we undertake the resummed computation of the gapped fraction 1Z which should cure the problem 
with the fixed-order computation we have just seen. Our calculation is based on the approach by Banfi-Marchcsini- 
Smye (BMS). We first describe their approach and adapt it for the problem at hand, and then present the numerical 
results. 



A. The BMS equation 

Consider a pair of cones pointing back-to-back in the direction of the beam axis (z-axis) as in Fig. |2j The opening 
angles of the right cone 8r and the left cone 0l need not be the same. A quark and an antiquark in the color singlet 
state ('dipole') are contained in the cones and moving in the directions f2 = (0,<fi) and il' = (6',<f>'), respectively, 
with the transverse momentum px- Let P t (0r, 9l, f2, fi') be the probability that the total amount of energy, or the 
transverse momentum, emitted from the dipole into the region outside the cones is less than E out - The 'evolution 
parameter' t is related to E out - Taking into account the running of the coupling, we havtj^] [151155] : 

f PT dk T a s (k T )N c 1 fa s {E out )\ 

where a s (pr) = %/(2bN c Iii(pt/Aqcd)), and b = (11N C — 2n/)/12JV c with rif = 5 being the number of active flavors. 

To zeroth order, P = 1. Within the large- iV c approximation, BMS derived a nonlinear integro-differential equation 
in t [16 which resums all the single logarithms (both Sudakov and non-global) (a s hipT/E ou t) ~ t™ which appear in 
the weak coupling expansion of P. The numerical solution to this equation has been previously obtained in |161 117] . 
Here we solve it again for the kinematic range relevant to the ATLAS data. In doing so, we switch to the rapidity 
variable 

Q 

y = In cot-, (4) 



5 With a fixed coupling approximation, we find r = (a a N c /n) \n(px/E ou t). 
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Figure 2. Br/l (vr/l = In cot 8 R /l/2) is the opening angle (rapidity) of the right/left cone. 



instead of the polar angle 9 in the arguments of P. Also, we suppress the dependence on the azimuthal angle <j) because 
P is a function of the difference A0 = — and below we need solutions only at A<f> = tt. Since the probability P 
must not depend on a specific Lorentz frame, we can evaluate it in the di-jet c.m.s. 

D / i\ o (Vr-Vl Vr~ Vl Vr + Vl , Vr + Vl\ ftC . 
P T {yR,yL,y,y) = Prl — ^ — , ^ — ,y ^ — ,y ^ — I- (5) 

Eq. ([5]) allows us to restrict ourselves to the symmetric configuration of the cones, which will greatly facilitate the 
computations below. 



B. The jet-veto cross section in hadron-hadron collisions 

As already suggested by BMS, one can utilize the above probability P T to constrain the inter-jet radiation in hadron- 
hadron collisions. The initial and final state soft radiations from a given partonic subprocesses may be deconstructed, 
in the large--ZV c approximation, into those from elementary dipoles by inspecting the flow of color. Consider the 
simplest P1P2 —> P3P4 process: qq' — > qq' or qq' — > qq' (different quark flavors). The partonic cross section is given by 

dt lD7rs z 

where s, t, and it are the standard partonic Mandclstam variables (s = X1X2S, etc) and (C F = (N? - l)/2N c = |) 

At large- N e , and in the one-gluon exchange, color flows as 1 — > 4 and 2 — >■ 3. Thus the radiation from the four-parton 
system (1234) factorizes into that from two dipoles (14) and (23). We require that the amount of energy emitted in 
the central region bounded by the edges of the jets p^ and p4 with jet radius R is less than E out (see, Fig. [2|. The 
cross section with this requirement is given by [16 



J veto 
a(J qq> 



(h A (§,i, u)P T (y3 - R,y4 + R,oo,y 4 )P T (y 3 - R,y 4 + R,y 3 , -oo) 



dt 167TS 2 

+h A (s,u,t)P T (y 3 - R,y 4 + R,oo,y 3 )P T (y 3 — R,y 4 + R,y 4 , -oo) J . (8) 



Note that we have added the 'it-channel' term although there is no it-channel diagram for the process qq' — > qq', 
because, when calculating the gapped cross section, we must sum over the two cases Ay = y 3 — y 4 > and Ay < 
since q and q' are not identical particles. Adding these two contributions effectively amounts to adding the it-channel 
term but restricting to the Ay > case, facilitating later computation. 
By boost invariance we have the relations 

PAys -R,y 4 + R, oo, y 4 ) = P T - R, + R, oo, -^j = P u {t, Ay) , 

P T (y 3 -R,y4 + R,y 3 ,~^) = Pr (^--R,-^ + R,^,-oc\ =P 23 (r,Ay), (9) 
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Figure 3. Comparison of the resummed veto fraction with the ATLAS measurement, for a fixed veto energy of E out = 20 GeV, 
in different bins of pr- The inner (green) uncertainty band is obtained taking into account only the renormalization and 
factorization scale uncertainties, while the outer (yellow) band also includes the subleading logarithmic uncertainty. For the 
ATLAS data, circles represent the case where the two leading jets are selected while the one where the most forward and 
backward jets are selected are represented by crosses. 



P T (y 3 -R,y 4 + R,oo,y 3 )=P T - R, + R, oo, ^ = P 13 { T ,Ay), 

P T (y3~R,y4 + R,y4,-™) = Pr - R,-~ + R,-~,-ocj =P 24 (r,A y ). (10) 
Clearly, Pi 4 = P23 and Pi 3 = P24, so that ^ takes the simpler form 

— TP" = T7T^2 ( hA ( S > t «) p i4^2 3 + h A (s, u, t)P 13 P 2i ) 

= ^ (h A (§, i, u)Pl + h A (s, u, f)i*) . (11) 

In Appendix we list the gapped cross section of all the other partonic subprocesses obtained in a similar way. The 
results are convoluted with the PDFs to give the hadronic gapped cross section 

V / dYxifi(x 1 ,p T )x2fj(x2,p T ) , (12) 

Jy min T at 



where Y = and ^max = -^min = 4.4 for the ATLAS setup. 



C. Results 



We compute the resummed gapped fraction 1Z by dividing ( 12 ) by the leading-order cross section (that is, the same 
expression except that P's are set to 1 everywhere). We have kept n/ = 5 as adequate in the kinematic range under 
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Figure 4. Comparison of the resummed veto fraction with the ATLAS measurement, for a fixed veto energy of E out = 20 GeV, 
in different bins of pr- See Fig. [3] for details. 

consideration and fixed the running of the coupling by imposing a s {Mz) — 0.12. The theoretical uncertainty on our 
predictions comes from different sources: First, the renormalization and factorization scale uncertainties are obtained 
by varying the scale of a s in the matrix elements and in the definition of r (Eq. ([3])) — in the latter case, we vary 
directly the scale of a s (kx) in the integrand — as well as the scale pt in the PDFs by a factor of 2 up and down. 
We exclude the cases where the renormalization scale /i r and the factorization scale \if are (fi r , fif) = (px/2, 2px) or 
(2pr,pr/2). It is interesting to note that the factor a^(/i r ) appearing in the matrix elements cancels when taking 
the cross-section ratio and thus the only renormalization scale uncertainty that remains is the one in the definition 
of r. Then, the uncertainty related to subleading logarithmic corrections that are not resummed by the leading 
a™ \n n (pt / E out ) series is estimated by varying the upper integration bound in the definition of r, Eq. ^ by a factor 
of 2 up and down. Finally, the uncertainty band is then taken as the envelope of all the resulting curves. 

The results are compared against the ATLAS measurements as published in [1 . Since we use the LO cross section 
where there are only two jets in the final state, we cannot distinguish the two ATLAS data sets based on different 
di-jet selection criteria. While naively we expect our predictions to be in better agreement with the "two-leading jets" 
selection, we shall see below that we better reproduce the data based on the most forward/backward jets. 

Fig. [3] shows the dependence on the jet rapidity separation for different bins in jet px- We note that the model is 
overall in good agreement with the measurement. Needless to say, the unphysical behavior, 1Z > 1 or 1Z < 0, of the 
fixed-order computation (see Fig. [I]) has disappeared. In Fig. [4] the result is presented over the jet pr for different bins 
in Ay. In both cases, the agreement with the ATLAS data is good, especially with the sample obtained by selecting 
the most forward/ backward jets. The predictions systematically undershoot the data based on the two leading jets 
although they remain in agreement once the uncertainty due to subleading effects are taking into account. 

Then, Fig. [5] demonstrates the dependence on the jet veto threshold energy E out . We note first that in our results 
the veto fraction saturates to unity as the threshold -Eout approaches the jet px scale, as expected. In the data, it 
falls short of unity because of the NLO (2 — > 3) corrections. On the other hand, the agreement is best for the smaller 
values of i? ou t, which is expected as well since the formalism we use requires a large scale separation between the jet 
Pt and the threshold -E ut- 

Finally, it is worth mentioning that uncertainties in the choice of the PDF set mostly cancel in the ratio ([lj. We 
checked this by trying two different sets, MRST2002 [27J and CT10 [25] , and found no noticeable difference. 

IV. CONCLUSIONS AND OUTLOOK 

In this paper, we have investigated the QCD resummation of the Sudakov and the non-global logarithms induced 
by soft gluon emissions in the context of the jet veto cross-section. We conclude that the ATLAS measurement 
is well described by tree-level QCD supplemented with the jet veto probability calculated perturbatively using the 
BMS equation which resums both logarithms mentioned above. Actually, the impact of the non-global logarithms is 
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Figure 5. Comparison of the resummed veto fraction with the ATLAS measurement, for different kinematic bins, as a function 
of the veto threshold E ou t. See Fig. [3] for details. 



modest compared to that of the Sudakov logarithms: with E out as large as 20 GeV, we estimate that the non-global 
contribution reduces the gap fraction by about 15%. 

Having said this, we must comment on uncertainties in our results other than the scale uncertainties already 
examined. Above all, our description does not allow to disentangle between the two methods proposed by ATLAS 
to identify the di-jet system. While our approach better reproduces the data where the most forward and the most 
backward jets are selected, it is quite below the measurement obtained by selecting the two hardest jets in the event. 
This somewhat goes against the naive expectation that our description should better reproduce the latter situation. 
This probably means that effects other than the ones included in the BMS equation are at play. Still, we have shown 
that the resummation of the soft gluon emissions gives a reasonable description of the jet veto probability. 

Various additional effects would play a role if we wanted to improve our predictions. First, our approach resums 
the soft gluon emissions at the leading logarithmic accuracy, so subleading effects would potentially be important. 
[A rough estimate of this is shown by the yellow band in the figures.] Then, we do not include 1/N C corrections 
neither in the Sudakov logarithms nor in the non-global ones. The 1/N C corrections in the Sudakov logarithms are 
fully taken into account in [TH], but those for the non-global ones pose a serious theoretical challenge. Also, the NLO 
(2 — > 3) corrections to the hard parton cross sections can be important, especially at small Ay and for discriminating 
the two ATLAS data sets corresponding to two different definitions of the di-jet system. In principle, our resummed 
approach can be extended to 2 — > 3 processes, but the flow of color will be considerably more complicated and it is 
not clear to us at the moment if such an extension is practical. Alternatively, it may be more practical to match the 
BMS predictions to the fixed-order predictions from Section II which are expected to provide more accurate NLO 
corrections at least in the small Ay region. Finally, we have seen that the data are well described by the color octet 
contribution alone, without introducing additional BFKL-like (singlet) contributions. As already mentioned in the 
introduction, the ATLAS choice E out ^S> ^qcd significantly reduces the sensitivity to BFKL exchanges. Presumably, 
we need to go to lower values of E out [H], thus approaching the case of a perfect gap (E out — > 0) [HEIGH], to enter 
the BFKL-dominated regime. It is nevertheless interesting to mention that these corrections would tend to increase 
our results for the gap fraction TZ, especially at large Ay: Our estimates for 1Z based on a pure-BFKL calculation are 
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in the range 0.8-0.9, and the HEJ event generator, which includes BFKL effects, also predicts values larger than the 
ones obtained by ATLAS. A combination of the BMS and BFKL contributions could thus bring our predictions in 
better agreement with the ATLAS measurement obtained based on the two hardest jets in the event. We leave this 
for future study. 
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Appendix A: Parton subprocess 

In this Appendix we list the gapped cross section of all the partonic subprocesses which contribute to the sum in 



(12 1. Some of them were already considered in |16j . The others are simply obtained by crossing symmetry. We shall 



need the following building blocks of the cross section [TH1 1ST] 



h B (s,t,u) = g i — — _ + _- , (Al) 



N c V t 2 N c t 

^( S ,M)=^|(* 2 ^ 2 -i), (A2) 

, n / \ „ a ( tu su u 2 \ , . . 

^(M,.) = 2. 9 ^^l---- + -j, (A3) 

and the solutions to the BMS equation for different dipole configurations 

Pr(y 3 -R,yA + R,y 3 ,yi)=P T - R, ~ + R ^~^) = P3a(t,A v ), (A4) 

P T {y s - R,yi + R,oc,-oo) = P T - R,-^- + R,oo, -oo) = P 12 (r,Ay). (A5) 
Note that P 34 is evaluated at Acj) = 4> 3 — 04 = n as appropriate for back-to-back jets. 

• <Z<f QQ' (different quark flavors) 
This can be obtained from qq' — > qq' by crossing s ^ u and 2 -H- 4 



J veto 
a(J qq> 



di 16ns 



l — (h A (u, t, s)P 12 P 34 + h A (t, u, s)Pi 2 P 34 ) , (A6) 



where we added a w-channel diagram (by the same reason as in ([8|) which is obtained by crossing t <H> u, 3 «4 4. 



• qq q'q' (different quark flavors) 

Obtained from qq' — > qq' via crossing s H t, 2 O 3 



di 16tts 2 



^ (h A {t, s, u)P u P 23 + h A (u, s, t)P 13 P 2 . 



1&TTS 2 

qq — > qq, qq — > qq (the same flavor) [IB] : 



' 'h A (t, s, u)Pl + h A (u, s, t)P 2 3 ) . (A7) 



d(7 qq 



: — K^(h B (s 7 i,u)P 14 P 23 + h B (s,u,i)P 13 P 24 ) 
dt lbns z V / 

= ") p i 2 4 + hB ( § > fi. fyPxn) ■ (A8) 
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Compare with (11). Since the final state particles are identical, one integrates over only half of the phase 
space, that is, only the Ay > case. Note that the interference term between the two diagrams for qq — > qq is 
subleading in N c and does n ot co rrespond to a process with definite color flow. It was distributed symmetrically 



into the second term of h B (All [3T]. Similar comments apply to the other subleading contributions in h and 
hP below. 

• qq -> qq- 

Obtained from qq — >• qq by crossing s O u, 2 O 4 and adding the Ay < contribution (f f} «, 3 fl 4) because 
of nonidcntical particles in the final state 



j veto -t 
qq-^qq ' 



h B (u, t, s)P 12 P 34 + h B (u, a, t)P 13 P 24 + h B (t, u, s)P 12 P 34 + h B (t, s, u)P 14 P 2 . 



di 167TS 2 

1 A2 ( hB (^ *> s)P 12 P 34 + h B (u, 8, t)Pi 2 3 + h B {l u, s)P 12 P 34 + h B (t, 8, u)P? 4 ) . (A9) 
loirs- 2 V / 



• qq-+ gg 



J veto 
u,u 

"23 



dt 167TS 

1 

167TS 



l — (h C ( S , t, u)P 34 P 13 P 24 + h C ( S , U, t)P 34 P 14 P 2 , 
ITS'* \ 

^2 (/i c (s, t, u)P 34 Pl 3 + h c (s, u, i)P 34 P? 4 ) . (A10) 



99 qq- 

The same as qq — > gg except for the color factor (3/8) 2 and the addition of the Ay < contribution (t «-» u, 
3 o 4) 



J veto 1 /n\2 

di 16tts 2 



(/i c (s, t, u)P 34 Pl 3 + h c (s, u, i)P 34 P? 4 ) x 2 . (All) 



• qq q9 an d qg ~ ^ w- 

Obtained from qq — > gg by crossing s t, 2 <h> 3, multiplying the color factor 3/8, adding the Ay < 
contribution (t O it, 3 <-» 4) 

rfa veto —13/ \ 
= T7T^2 o [ hC (i> § >u)P 24 P 12 P 34 + h c (t,u,s)P 24 P 14 P 23 + h c (u,s,t)P 23 P 12 P 34 + h c (u,t,s)P 23 P 13 P 24 ) 
dt lons z 8 V / 

= T7T^2 ^ ( hC & s > u)P 13 P 12 P U + h c (t, u, s)P 13 Pf 4 + h c (u, s, t)P 14 P 12 P 34 + h c (u, t, s)P 14 P? 3 ) . (A12) 

1D7TS O V / 

Note the overall minus sign. 

• 99^ 99 US]: 

da veto 1 / 

9 "7" = T7T^2 *> u)P 12 P 13 P 24 P 34 + h D (§, u, t)P 12 P 14 P 23 P 34 + h D (u, t, s)P 14 P 24 P 13 P 23 

dt Ibns'' V 

= T7T^2 ( h °^ *> u)P 12 Pl 3 P 34 + h D (s, u, t)P 12 Pl 4 P 34 + h D {u, t, §)Pl \Pl 3 ) . (A13) 



[1] G. Aad et al. [ATLAS Collaboration], JHEP 1109, 053 (2011). 

[2] B. Abbott et al. [DO Collaboration], Phys. Lett. B 440, 189 (1998). 

[3] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 80, 1156 (1998). 

[4] F. Chevallier, O. Kepka, C. Marquet and C. Royon, Phys. Rev. D 79, 094019 (2009). 

[5] O. Kepka, C. Marquet and C. Royon, Phys. Rev. D 83, 034036 (2011). 



10 



[6] M. Bahr, S. Gieseke, M. A. Gigg, D. Grellscheid, K. Hamilton, O. Latunde-Dada, S. Platzer and P. Richardson et al, Eur. 
Phys. J. C 58 (2008) 639. 

T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605 (2006) 026. 

M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307 (2003) 001. 
S. Alioli, P. Nason, C. Oleari and E. Re, JHEP 1006 (2010) 043. 

J. R. Andersen and J. M. Smillie, Nucl. Phys. Proc. Suppl. 205-206 (2010) 205; JHEP 1106 (2011) 010. 
G. Oderda and G. F. Sterman, Phys. Rev. Lett. 81, 3591 (1998). 
J. Forshaw, J. Keates and S. Marzani, JHEP 0907, 023 (2009). 
M. Dasgupta and G. P. Salam, Phys. Lett. B 512, 323 (2001). 
E. Avsar, Y. Hatta and T. Matsuo, JHEP 0906, 011 (2009). 
A. Banfi, G. Corcella and M. Dasgupta, JHEP 0703, 050 (2007). 
A. Banfi, G. Marchesini and G. Smye, JHEP 0208, 006 (2002). 
Y. Hatta and T. Ueda, Phys. Rev. D 80, 074018 (2009). 

R. M. Duran Delgado, J. R. Forshaw, S. Marza ni and M. H. S eymour, JHEP 1108, 157 (2011). 
M. Deak, F. Hautmann, H. Jung and K. Kutak, [arXiv:1112.638"6l [hep-ph], 
M. Deak, F. Hautmann, H. Jung and K. Kutak, Eur. Phys. J. C 72, 1982 (2012). 
S. Alioli, J. R. Andersen, C. Oleari, E. Re and J. M. Smillie, Phys. Rev. D 85, 114034 (2012). 
E. Gerwick, T. Plehn, S. Schumann and P. Schichtel, JHEP 1210, 162 (2012). 
M. Cacciari, G. P. Salam and G. Soyez, JHEP 0804 (2008) 063. 
G. Soyez, Phys. Lett. B 698, 59 (2011). 

Z. Nagy, Phys. Rev. Lett. 88 (2002) 122003; Phys. Rev. D 68 (2003) 094002. 
M. Cacciari and G. P. Salam, Phys. Lett. B 641 (2006) 57; 
M. Cacciari, G. P. Salam and G. Soyez, Eur. Phys. J. C 72 (2012) 1896. 

A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28 (2003) 455. 

G. Marchesini and A. H. Mueller, Phys. Lett. B 575, 37 (2003). 

H. -L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsk y, J. Pumplin and C. -P. Yuan, Phys. Rev. D 82 (2010) 074024. 
C. Marquet, C. Royon, M. Trzebinski and R. Zlebcik, [arXiv:1212.205"9l [hep-ph], 
G. Marchesini and B. R. Webber, Nucl. Phys. B 310, 461 (1988). 



